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Abstract. In this paper we describe a singly exponential algorithm for com- 
puting the first Betti number of a given semi-algebraic set. Singly exponential 
algorithms for computing the zero-th Betti number, and the Euler-Poincare 
characteristic, were known before. No singly exponential algorithm was known 
for computing any of the individual Betti numbers other than the zero-th 
one. We also give algorithms for obtaining semi-algebraic descriptions of the 
semi-algebraically connected components of any given real algebraic or semi- 
algebraic set in single-exponential time improving on previous results. 



1. Introduction 

Let R be a real closed field and 5* C R'" a semi-algebraic set defined by a 
quantifier-free Boolean formula with atoms of the form P>0, P<0, P = 
for P € V C R[Xi, . . . , Xk]- We call S a T'-semi-algebraic set. If, instead, the 
Boolean formula has atoms of the form P = 0, P > 0, P < 0, P € V, and addition- 
ally contains no negation, then we will call S a P-closed semi-algebraic set. It is 
well known [18, 19, 17, 21, 1, 12] that the topological complexity of S (measured 
by the various Betti numbers of S) is bounded by 0{sd)'', where s = #(P) and 
d = maxpg-p deg(P). More precise bounds on the individual Betti numbers of S ap- 
pear in [2]. Even though the Betti numbers of 5* are bounded singly exponentially 
in fc, there is no singly exponential algorithm for computing the Betti numbers of S. 
This absence is related to the fact that there is no known algorithm for producing 
a singly exponential sized triangulation of S (which would immediately imply a 
singly exponential algorithm for computing the Betti numbers of S). In fact, the 
existence of a singly exponential sized triangulation, is considered to be a major 
open question in algorithmic real algebraic geometry. Moreover, determining the 
exact complexity of computing the Betti numbers of semi-algebraic sets is an area of 
active research in computational complexity theory, for instance counting versions 
of complexity classes in the Blum-Shub-Smale model of computation (see [9]). 

Doubly exponential algorithms (with complexity (sd)^°* ' ) for computing all 
the Betti numbers are known, since it is possible to obtain a triangulation of S 
in doubly exponential time using cylindrical algebraic decomposition [11, 6]. In 
the absence of a singly exponential time algorithm for computing triangulations of 
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semi-algebraic sets, algorithms with single exponential complexity are known only 
for the problems of testing emptiness [20, 4], computing the zero-tli Bctti number 
(i.e. the number of semi-algebraically connected components of S) [14, 10, 13, 5], 
as well as the Euler-Poincarc characteristic of S [1] . 

In this paper we describe the first singly exponential algorithm for computing 
the first Betti number of a given semi- algebraic set S C R''. In the process, we 
also give efficient algorithms for obtaining semi-algebraic descriptions of the semi- 
algebraically connected components of a given real algebraic or semi-algebraic set. 
These algorithms have complexity bounds which improve the complexity of the best 
previously known algorithm [16]. 

The rest of the paper is organized as follows. 

There are several ideas involved in the design of our algorithm for computing 
the first Betti number of a given semi- algebraic set, which corresponds to the main 
steps in our algorithm. We describe each of them separately in different sections. 

In Section 3 we recall the notion of a roadmap of an algebraic set [5] and indicate 
how to use it to construct connecting paths in basic semi-algebraic sets. 

In Section 4 we define certain semi-algebraic sets which we call parametrized 
paths and prove that under a certain hypothesis these sets arc semi-algebraically 
contractible. We also outline the input, output, and complexity of an algorithm 
computing a covering of a given basic semi-algebraic set, S C R*", by a singly 
exponential number of parametrized paths. 

In Section 5, we use the properties of parametrized paths proved in Section 4 
to give an algorithm (Algorithm 2) for computing a covering of a given closed and 
bounded semi-algebraic set by a single exponential sized family of closed, bounded 
as well as contractible semi-algebraic sets. The complexity of this algorithm is 
singly exponential. 

In Section 6, we recall some results from algebraic topology which allows us to 
compute the first Betti number of a closed and bounded semi-algebraic set from 
a covering of the given set consisting of closed, bounded and contractible sets. 
The main tool here is a spectral sequence associated to the Mayer- Vietoris double 
complex. We show how to compute the first Betti number once we have computed 
a covering by closed contractible sets and the number of connected components of 
their pair-wise and triple-wise intersections of the sets in this covering and their 
incidences. If the size of the covering is singly exponential, this yields a singly 
exponential algorithm for computing the first Betti number. Extensions of these 
ideas for computing a fixed number of higher Betti numbers in singly exponential 
time is possible and is reported on in a subsequent paper [3]. In Section 7, we 
describe an algorithm for computing the first Betti number of a given P-closed 
semi-algebraic set. 

In Section 8 we recall a technique introduced by Gabrielov and Vorobjov [12], for 
replacing any given semi-algebraic set by one which is closed and bounded and has 
the same homotopy type as the given set. In fact, we prove a slight strengthening 
of the main result in [12], in that we prove that the new set has the same homotopy 
type as the given one, while the corresponding result (Lemma 5) in [12] states that 
only the sum of the Betti numbers is preserved. The above construction allows 
us to reduce the case of general semi-algebraic sets to ones which are closed and 
bounded treated in Section 6 without any significant worsening of complexity. 



COMPUTING THE FIRST BETTI NUMBER 



3 



In Section 9, we describe an algorithm for computing the first Betti number of 
a general semi-algebraic set, using the construction described in Section 8 to first 
reduce the problem to the 7^-closed case already treated in Section 7. 

Finally, in Section 10 we indicate that the algorithms described in Section 4 
actually produces descriptions of the connected components of a given algebraic or 
semi-algebraic set in an efficient manner. 

2. Preliminaries 
Let R be a real closed field. For an element a e R wc let 

{0 if a = 0, 
1 if a > 0, 
-1 ifa<0. 

If P is a finite subset of R[Xi, . . . , Xk], we write the set of zeros of V in R*^ as 

Z{r,R'') ^{xeR'' \ /\ P{x) = 0}. 

Pev 

We denote by i?(0, r) the open ball with center and radius r. 

Let Q and V be finite subsets of R[Xi, ...,Xk], Z = Z(Q, R''), and Zr ^ Z r\ 
B{0,r). A sign condition on V is an element of {0, 1, —1}^. The realization of the 
sign condition a over Z, TZia, Z), is the basic semi-algebraic set 

{x e RM A ^(^) = ^ A sign(P(x)) = a{P)}. 

Q£Q Pev 

The realization of the sign condition a over Zj., Zr), is the basic semi-algebraic 
set TZ{a,Z) n B{0,r). For the rest of the paper, we fix an open ball B{0,r) with 
center and radius r big enough so that, for every sign condition a, TZ^a, Z) and 
TZ{<J, Zr) are homeomorphic. This is always possible by the local conical structure 
at infinity of semi-algebraic sets ([7], page 225). 

A closed and bounded semi-algebraic set 5* C R*^ is semi-algebraically triangula- 
ble (see [6]), and we denote by Hi{S) the i-th simplicial homology group of S with 
rational coefficients. The groups Hi{S) are invariant under semi-algebraic home- 
omorphisms and coincide with the corresponding singular homology groups when 
R = M. We denote by bi{S) the i-th. Betti number of S (that is, the dimension 
of Hi{S) as a vector space), and by b{S) the sum J2i bi{S). For a closed but not 
necessarily bounded semi-algebraic set S" C R , we will denote by Hi{S) the i-th 
simplicial homology group of fl B{Q,r), where r is sufficiently large. The sets 
S n -6(0, r) are semi-algebraically homeomorphic for all sufficiently large r > 0, by 
the local conical structure at infinity of semi-algebraic sets, and hence this definition 
makes sense. 

The definition of homology groups of arbitrary semi-algebraic sets in R'" requires 
some care and several possibilities exist. In this paper, we define the homology 
groups of realizations of sign conditions as follows. 

Let R denote a real closed field and R' a real closed field containing R. Given 
a semi-algebraic set in R*^, the extension of S to R', denoted Ext(S', R'), is the 
semi-algebraic subset of R defined by the same quantifier free formula that defines 
S. The set Ext(5', R') is well defined (i.e. it only depends on the set S and not on 



4 



SAUGATA BASU, RICHARD POLLACK, AND MARIE-FRANgOISE ROY 



the quantifier free formula chosen to describe it). This is an easy consequence of 
the transfer principle [6]. 

Now, let S C R*" be a T'-semi-algebraic set, where V = {Pi, . . . , Ps} is a finite 
subset of R[Xi, . . . , Xk]- Let <f>{X) be a quantifier- free formula defining S. Let 
Pi = where the a^.Q G R. Let A ~ {. . . , Ai^a, ■ ■ ■) denote the vector 

of variables corresponding to the coefficients of the polynomials in the family V, 
and let a = (. . . , a^^a, . . .) S R^ denote the vector of the actual coefficients of 
the polynomials in V. Let tp{A, X) denote the formula obtained from (f>(X) by 
replacing each coefficient of each polynomial in V by the corresponding variable, so 
that ipiX) = ■i/;(a,X). It follows from Hardt's triviality theorem for scmi-algebraic 
mappings [15], that there exists, a' £ R^g such that denoting by S' C M^ig the 
semi-algebraic set defined by ip{a',X), the semi- algebraic set Ext(S",R), has the 
same homeomorphism type as 5*. We define the homology groups of S to be the sin- 
gular homology groups of Ext(S", R). It follows from the Tarski-Seidenberg transfer 
principle, and the corresponding property of singular homology groups, that the 
homology groups defined this way are invariant under semi-algebraic homotopies. 
It is also clear that this definition is compatible with the simplicial homology for 
closed, bounded scmi-algcbraic sets, and the singular homology groups when the 
ground field is R. Finally it is also clear that, the Betti numbers are not changed 
after extension: bi{S) = 6j(Ext(S', R')). 

3. ROADMAP OF A SEMI- ALGEBRAIC SET 

We first define a roadmap of a scmi-algcbraic set. Roadmaps are crucial in- 
gredients in all singly exponential algorithms known for computing connectivity 
properties of semi-algebraic sets such as computing the number of connected com- 
ponents, as well as testing whether two points of a given semi-algebraic set belong 
to the same semi-algebraically connected component. 

Wc use the following notations. Given x = (xi , . . . , x^) we write Xi for (xi , . . . , Xi), 
and Xi for (x^+i, . . . We also denote by TTi...j the projection, x i~> Xj. Given a 
set 5" C R'', 2/ e R-' wc denote hy Sy = S f] iT^^,j{y). 

Let S* C R*^ be a semi-algebraic set. A roadmap for 5 is a semi-algebraic set 
M of dimension at most one contained in S which satisfies the following roadmap 
conditions: 

• RMi For every semi-algebraically connected component of 5*, D fl M is 
semi-algebraically connected. 

• RM2 For every a; S R and for every semi-algebraically connected component 
D' oiS^, D'r\M^%. 

Wc describe the construction of a roadmap M for a bounded algebraic set 
Z{Q,Y&) which contains a finite set of points N oi'L{Q,Y&). A precise description 
of how the construction can be performed algorithmically can be found in [6] . 

A key ingredient of the roadmap is the construction of a particular finite set 
of points having the property that, they intersect every connected component of 
Z((5, R*^). We call them Xi-pseudo-critical points, since they are obtained as limits 
of the critical points of the projection to the Xi coordinate of a bounded nonsingu- 
lar algebraic hypcrsurfacc defined by a particular infinitesimal deformation of the 
polynomial Q. Their projections on the Xi-axis are called pseudo-critical values. 
These points are obtained as follows. 
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We denote by R(C) the real closed field of algebraic Puiseux series in C, with 
coefHcients in R [6]. The sign of a Puiseux series in R(C) agrees with the sign of 
the coefficient of the lowest degree term in ^. This induces a unique order on R(C) 
which makes ( infinitesimal: ^ is positive and smaller than any positive clement of 
R. When a e R(C) is bounded by an element of R, lim^(a) is the constant term of 
a, obtained by substituting for ^ in a. We now define the deformation Q of Q as 
follows. Suppose that Z{Q,R'') is contained in the ball of center and radius 1/c. 
Let d be an even integer bigger than the degree d oi Q, 

(3.1) Gk{d, c) = c^ixf +--- + Xi+X^ + --- + Xi)- {2k - 1), 

(3.2) Q^CGkid,c) + {l-C)Q. 

The algebraic set Z(Q,R(C) ) is a bounded and non-singular hypersurface lying 
infinitesimally close to Z((5, R ), and the critical points of the projection map onto 
the Xi co-ordinate restricted to Z{Q, R(C) ) form a finite set of points. We take the 
images of these points under lim^ and we call the points obtained in this manner 
the Xi-pseudo-critical points of Z((5,R'°). Their projections on the Xi-axis are 
called pseudo-critical values. 

The construction of the roadmap of an algebraic set containing a finite number 
of input points J\f of this algebraic set is as follows. We first construct X2-pseudo- 
critical points on Z(Q, R'°) in a parametric way along the Xi-axis, by following con- 
tinuously, as X varies on the Xi-axis, the X2-pseudo-critical points on Z{Q,K )x- 
This results in curve segments and their endpoints on Z((5, R ). The curve segments 
arc continuous semi-algebraic curves parametrized by open intervals on the Xi-axis, 
and their endpoints are points of Z((5,R'') above the corresponding endpoints of 
the open intervals. Since these curves and their endpoints include, for every a; G R, 
the X2— pseudo-critical points of Z((5, R'^)x! they meet every connected component 
of Z{Q,R'^)x- Thus the set of curve segments and their endpoints already satisfy 
RM2. However, it is clear that this set might not be semi-algebraically connected 
in a semi-algebraically connected component, so RMi might not be satisfied. We 
add additional curve segments to ensure connectedness by recursing in certain dis- 
tinguished hyperplanes defined hy Xi = z for distinguished values z. 

The set of distinguished values is the union of the Xi-pseudo-critical values, the 
first coordinates of the input points Af and the first coordinates of the endpoints 
of the curve segments. A distinguished hyperplane is an hyperplane defined by 
Xi ~ V, where w is a distinguished value. The input points, the endpoints of the 
curve segments and the intersections of the curve segments with the distinguished 
hyperplanes define the set of distinguished points. . 

Let the distinguished values he vi < . . . < vg. Note that amongst these are the 
Xi-pseudo-critical values. Above each interval (vi, Wi+i), we have constructed a col- 
lection of curve segments Ci meeting every semi-algebraically connected component 
of Z(Q,R'=)^, for every v € {vi,v.i+i). Above each distinguished value Vi, we have 
a set of distinguished points Afi . Each curve segment in Ci has an endpoint in A/i 
and another in A/i+i. Moreover, the union of the A/i contains M. 

We then repeat this construction in each distinguished hyperplane Hi defined by 
Xi = Vi with input Q{vi, X2, . . . , Xk) and the distinguished points in A/i. Thus, we 
construct distinguished values, Vi^i, . . . , Wi^£(i) of Z{Q{vi, X2, ■ ■ ■ , Xk), R (with 
the role of Xi being now played by X2) and the process is iterated until for / = 
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(ii, . . . ,ifc_2), 1 < ii < ^, < ik-2 < ■ ■ ■ ,ik~3), we have distinguished 

values W74 < . . . < vij(^j-) along the Xk^i axis with corresponding sets of curve 
segments and sets of distinguished points with the required incidences between 
them. 

The following proposition is proved in [5] (see also [6]). 

Proposition 1. The semi-algebraic set RM(Z((5, R'''),7V) obtained by this construc- 
tion is a roadmap for Z{Q, R'^') containing J\f. 

Note that if a; e Z(Q, R*"). RM(Z((5, R'^), {x}) contains a path, j{x), connecting 
a distinguished point p of RM(Z((3, R''")) to x. 

Later in this paper we shall examine the properties of parametrized paths which 
are the unions of connecting paths starting at a given p and ending at x, where 
X varies over a certain semi-algebraic subset of Z{Q,R''). In order to do so it is 
useful to have a better understanding of the structure of these connecting paths - 
especially, of their dependence on x. 

Recall that given x = {xi, . . . ,Xk) we write Xi for and Xi for 

(xi+i, . . . , Xk)- We first note that for any x = {xi, . . . , Xk) £ Z{Q,K^), we have 
by construction that, RM(Z((3, R*")) is contained in RM(Z((5, R*"), {x}). In fact, 

RM(Z(g, r'=), {x}) = RM(Z(Q, r'^)) U RM(Z(g, r''%, , M^J, 

where Mxi consists of xi and the finite set of points obtained by intersecting the 
curves in RM(Z((5, R'^)) parametrized by the Xi-coordinate, with the hyperplane 

^iH^i)- ^ ^ ' 

A connecting path 7(2;) (with non-self intersecting image) joining a distinguished 
point p of RM(Z(Q,R'')) to x can be extracted from RM(Z(Q, R''), {x}). The 
connecting path j{x) consists of two consecutive parts, 70(2;) and ri(x). The path 
7o(a;) is contained in RM(Z((5, R*^)) and the path ri(x) is contained in Z{Q,K'')x^- 
The part 70(2;) consists of a sequence of sub-paths, 7o,Oj ■ • ■ j 7o,m- Each 70, i is 
a semi-algebraic path parametrized by one of the co-ordinates Xi, . . . ,Xk, over 
some interval [ao,i, 60, i] with 70,0(00,0) = P- The semi-algebraic maps, 70,07 • ■ • j 70,™ 
and the end-points of their intervals of definition ao,o, ^o,Oj • ■ ■ i t^o,™! ^o,™ are all 
independent of x (upto the discrete choice of the path 7(2;) in RM(Z(Q, R*^), {a;})), 
except 60, m which depends on xi. 

Moreover, ri(a;) can again be decomposed into two parts, 71 (x) and T2{x) with 
r2(a;) contained in Z{Q,K'')x2 ^^'^ on. 

If y ~ (j/i, . . . , j/fc) e Z((5,R ) is another point such that xi 7^ yi, then since 
Z(Q,R'')j.j and Z{Q,K'^)y-^^ are disjoint, it is clear that 

RM(Z(Q, R'^), {x}) n RM(Z(Q, R'^), {y}) = RM(Z(Q, R^)). 

Now consider a connecting path "f{y) extracted from RM(Z((5, R*^), {y}). The im- 
ages of ri(a;) and ri(y) are disjoint. If the image of 70 (y) (which is contained in 
RM(Z(Q,R'=)) follows the same sequence of curve segments as 70(2;) starting at p 
(that is. it consists of the same curves segments 70, oj • ■ ■ j7o,rn as in 70(2;)), then it 
is clear that the images of the paths 7(2;) and 7(2/) has the property that they are 
identical upto a point and they are disjoint after it. 

More generally, if the points x and y are such that, Xi = y^,! < i < j and 
Xj^i ^ Uj+i, then the paths rj^i{x) and rj_|.i(y), contained in Z((3, R'^)^^^^ and 
Z((3, R'°)yj.^j^ respectively, will be disjoint. Moreover if the paths 70(2;), . . . ,7^(2;) 
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and jo{y),. . . ,7j(y) are composed of the same sequence of curve segments, then 
7(x) and 7(2/) wiU also have the divergence property. 

Wc now consider connecting paths in the scmi-algcbraic setting. We are given 
a polynomial Q e R[Xi, . . . ,Xk] such that Z(Q,R'^) is bounded and a finite set 
of polynomials V C D[Xi, . . . ,Xk] in strong ^-general position with respect to Q. 
This means that any £ + 1 polynomials belonging to V have no zeros in common 
with Q in R*", and any £ polynomials belonging to V have at most a finite number 
of zeros in common with Q in R*^ . 

For every point x of Z((5,R'^), we denote by a{x) the sign condition on V at x. 
Let 7^(CT(a;),Z(Q,R'=)) = {x £ Z{Q,R'') \ Ap^p sign(P(a;)) G CT(a;)(P)}, where a 
is the relaxation of a defined by 

= {0} if a = 0, 
W={Q,1} if (7 = 1, 
a = {0,-1} if cr = -l. 

Wc say that ^(a;) is the weak sign condition defined by x on V. Wc denote by 'P(x) 
the union of {Q} and the set of polynomials in V vanishing at x. 

The connecting algorithm associates to a; e Z(Q,R''') a path entirely contained 
in the realization of a{x) connecting a; to a distinguished point of the roadmap 
of some Z(7'',R'^), with V{x) C V. The connecting algorithm proceeds as follows: 
construct a path 7 connecting a distinguished point of RM(Z((5, R'^)) to x contained 
in RM(Z((3, R*^), {x}). If no polynomial of V \ V{x) vanishes on 7, we are done. 
Otherwise let y be the last point of 7 such that some polynomial of V\V{x) vanishes 
at y. Now keep the part of 7 connecting y to cc as end of the connecting path, and 
iterate the construction with noting that the realization of is contained in 
the realization of a{x)., and V \ V{y) is in £—1 strong general position with respect 
toZ(7^(y),R'=). 

As in the algebraic case, two such connecting paths which start with the same 
sequence of curve segments will have the divergence property. This follows from the 
divergence property in the algebraic case and the recursive definition of connecting 
paths. 

Formal descriptions and complexity analysis of the algorithms described above 
for computing roadmaps and connecting paths of algebraic and basic semi-algebraic 
sets can be found in [6] (Algorithm 15.12 and Algorithm 16.8). 

4. Parametrized paths 

We are given a polynomial Q G R[Ari, . . . ,Xk\ such that 7j{Q,T&) is bounded 
and a finite set of polynomials V G D[A"i, . . . , Xk\ in strong fc'-general position with 
respect to Q, where k' is the dimension oi7j{Q,T&). 

We show how to obtain a covering of a given 7^-closed scmi-algcbraic set con- 
tained in Z((5, R*^) by a family of semi-algebraically contractible subsets. The con- 
struction is based on a parametrized version of the connecting algorithm: we com- 
pute a family of polynomials such that for each realizable sign condition <j on this 
family, the description of the connecting paths of different points in the realization, 
7?,(tT, Z(Q, R*^)), are uniform. We first define parametrized paths. A parametrized 
path is a semi-algebraic set which is a union of semi-algebraic paths having the 
divergence property. 

More precisely, 
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Definition 4.1. A parametrized path 7 is a continuous semi-algebraic mapping 
from V C R*"'^^ R'', such that, denoting by U ^ Tri...k{V) C R*", there exists a 
semi-algebraic continuous function £ : U ^ [0, +00), and there exists a point a in 
R*", such that 

(1) V = {(.T,t) \ x€U,0<t< i{x)}, 

(2) VxeU, 7(3;, 0) = a, 

(3) Vx e U, jixj{x)) = X, 
(4) 

\/xe uyye uy se [0,£{x)],yte [OJ{y)] 
(7(2;, s) ='-f{y,t) ^s^t), 

(5) 

\f X e Uy y e Uy s e [0,mm{£{x),e{y))] 
(7(x, s) = -fiy, s) ^ V t < s j{x, t) = 7(2/, t)) . 

Given a parametrized path, 7 : F — > R*", we will refer to ?7 = 7ri...kiV) as its 
base. Also, any semi-algebraic subset U' C U of the base of such a parametrized 
path, defines in a natural way the restriction of 7 to the base U' , which is another 
parametrized path, obtained by restricting 7 to the set V C V, defined by V = 
{{x,t) \ xeU',0<t< e{x)}. 

Proposition 2. Let 7 : 1/ — > R'^ be a parametrized path such that U = T^i...k{V) is 
closed and bounded. Then, the image of 7 is semi-algebraically contractible. 

We thank A. Gabrielov and N. Vorobjov for pointing out an error in a previous 
version of this paper, where the same proposition was stated without any extra 
condition on U. In fact. Proposition 2 is not true if we do not assume that U is 
closed and bounded. 

Proof, (of Proposition 2) Let W = Im(7) and M = sup^g^; £{x). We prove that the 
semi-algebraic mapping 4> : W x [0, M] — > W sending 

• (7(x,t),s) to 7(x, s) if t > s, 

• (7(0;, t), s) to 7(x, t) if t < s. 

is continuous. Note that the map is well-defined, since 7(x, t) = 7(2;', t') ^ t ^ t', 
by condition (4). 
Since </> satisfies 

(j){-f{x,t),0) = a, 
(b{j{x,t),M)=^{x,t), 

this gives a semi-algebraic continuous contraction from W to {a}. 

Let w eW,s e [0, M]. Let e > be an infinitesimal, and let {w' , s') £ Ext(W x 
[0, M], R(e)) be such that limg{w' , s') = {w, s). In order to prove the continuity of 
(/) at w it suffices to prove that limg Ext(^, R{e)){w', s') = (p{w, s). 

Let w = ^{x,t) for some x G U,t G [0,£{x)], and similarly let w' = {x',t') for 
some x' e Ext(C/,R(e)) and t' G [0, Ext(£, R(e))(x')]. Note that lime(a;') G U since 
U is closed and bounded and lim^ t' G [0, €(lim£ x')]. 
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Now, 

^{x,t) = w 

= lime(i(;') 

= limeExt(7,R(e))(a;',f') 

= 7(limea;',limgt'). 

Condition (4) now implies that linij t' = t. 

Without loss of generality let t' > t. The other case is symmetric. We have the 
following two sub-cases. 

Case s' > t': Since s, i e R and lim^ s' = s and lim^ t' = t, wc must have that 
s > i. In this case Ext(0, R(e))(u;', s') = Ext(7, R(e))(.T', i')- Then, 

lim. Ext((/),R(e))(w',s') = lim. Ext(7, R(£))(x', i') 

= lim£ w' 

= w 

= (j)[w,s). 

Case s' < t': Again, since s, t G R and lim^ s' = s and lim^ t' = t, we must have 
that s < t. 

In this case we have, 

lime(^(w',s') = limsExt(7,R(£))(a;',s') 

— 7(lim£ a;', linig s') 

= 7(lim£ x', s). 

Now, 

7 (linig x' ,t) ~ 7 (limg x' , lim^ t' ) 

= limeExt(7,R(£))(a;',i') 

= lime w' 

~ w 

= l{x,t). 

Thus, by condition (5) we have that 7(lime x' , s") = 7(2;, s") for all s" < t. Since, 
s <t, this implies, 

limeExt(0,R(e))(w',s') = lim. Ext(7, R(£))(w', s') 

= 7(lime a;', lim^ s') 

= 7(x,s) 

= (l}{w,s). 

This proves the continuity of (j). □ 

We now describe how to compute parametrized paths in single exponential time 
using a parametrized version of the connecting algorithm. Wc describe the input, 
output and complexity of the algorithm which appears in [6] (Algorithm 16.15). 

Algorithm 1. (Parametrized Bounded Connecting) 
Input. 

— a polynomial Q e D[Xi, . . .,Xk], such that Z{Q,R'') C B{0, 1/c), 

— a finite set of polynomials V C r)[Xi, . . . , X^] in strong fc'-general 
position with respect to Q, where k' is the dimension of Z(Q,R'^). 

Output. 

— a finite set of polynomials A C R[Ari, . . . , X^], 
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— a finite set O of quantifier free formulas, with atoms of tlie form 
P = 0, P > 0, P < 0, P <E A, sucli tliat for every semi-algebraically 
connected component S of the reahzation of every weak sign condi- 
tion on V on Z(Q,R'''), there exists a subset 6(5) C 6 such that 

s= \J n{e,z{Q,R')), 

— for every 6 d Q, a, parametrized path 

with base Ue TZ{e, Z{Q, R'')), such that for each y e TZ{0, Z(Q, R'')), 
Im 7e(?/, ■) is a semi-algebraic path which connects the point y to a 
distinguished point ae of some roadmap RM{Z{T" U {Q},R'')) where 
V C "P, staying inside TZ(a{y),Z{Q,R'')). 
Complexity, s*"' +irfO(fc )^ where s is a bound on the number of elements 
of V and d is a bound on the degrees of Q and the elements of V. 

Proof of correctness. Given the proof of correctness of Algorithm 16.15 
(Parametrized Bounded Connecting) in [6], the only extra property that we need 
to prove is that for each S , 76 is a parametrized path. It is easy to see that 
satisfies the conditions of Definition 4.1, using the divergence property of the paths 
7(1/, ■) (see discussion in Section 3). □ 

5. Constructing coverings of closed semi-algebraic sets by closed 

contractible sets 

We are again given a polynomial Q G R[Xi,...,Xk] such that Z{Q,R'^) is 
bounded and a finite set of polynomials P C D[Xi, . . . , Xk] in strong fc'-general 
position with respect to Q, where k' is the dimension of Z{Q,R''). We describe an 
algorithm for computing closed contractible coverings of T'-closed semi-algebraic 
sets, using the results of Section 4. 

For the rest of this section we fix a T'-closed semi-algebraic set S contained 
in Z(Q,R'^) and let =f^A = t. We denote by Sign(yl, 5) the set of realizable sign 
conditions of A on Z(Q,R''') whose realizations are contained in S. We continue 
to follow the notations of Algorithm 1. For each a S Sign(^, S") TZ{<7,Z{Q,R )) is 
contained in 7?.(0, Z(Q, R*^)) for some 6 E O. We denote by 70- the restriction of 
7e to the base 7?.(cr, Z(Q, R'^))- Since Ti{a,Z{Q,R'^)) is not necessarily closed and 
bounded, Im 7o- might not be contractible. In order to ensure contractibility, we 
restrict the base of to a slightly smaller set which is closed, using infinitesimals. 

We introduce infinitesimals e2t ^ £21-1 ^ • • • ^ £2 ^ £1 > 0. For i = 1, . . . ,2t 
we will denote by R^ the field R(e2t) ■ • ■ (£i) and denote by R' the field Ri. 

For a e Sign(^, S) we define the level of a by, 

level(cr) = #{P e A I a{P) = 0}. 

Given a e Sign(^, S), with level((T) = j, we denote by TZ{(j^) the set defined on 
Z{Q,R2j) by the formula cr_ obtained by taking the conjunction of 

P = 0, for each P eA such that a{P) = 0, 
P > £2j, for each P £ A such that cr(P) = 1, 
P < — £2j, for each P G A such that a{P) ~ —1. 
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Notice that TZ{a-) is closed, bounded and contained in Ti{a, 'Zj{Q, R2j))- Proposi- 
tion 2 implies, 

Proposition 3. ^o-iT^i'^-)) is semi- algebraically contractible. 

Note that the sets 7cr(7?.((T_)) do not necessarily cover S. So we are going to 
enlarge them, preserving contractibility, to obtain a covering of S. 

Given a £ Sign(^, S), with leve\{a) = j, we denote by TZ{a^) the set defined on 
Z((3,R2j_i), by the formula obtained by taking the conjunction of 

— £2^-1 < P < £2j-i, for each P E A such that <t{P) — 0, 
P > £2j, for each P £ A such that cr(P) = 1, 
P < ~e2j, for each P G A such that cr(P) = — 1. 

with the formula (j) defining S. Let Co- be the set defined by, 

a =7a(7^(a_))U7^(a^). 

We now prove that 
Proposition 4. is semi-algebraically contractible. 

Let C be a closed and bounded semi-algebraic set contained in R(e)'^. We can 
suppose without loss of generality that C is defined over R[e] (see for example 
Proposition 2.80 in [6]). We denote by C* the semi-algebraic subset of R'^ defined 
by replacing e by t in the definition of C. Note that is nothing but C. 

We are going to use the following lemma. 

Lemma 5.1. Let B be a closed and bounded semi-algebraic set contained in I& and 
let C be a closed and bounded semi- algebraic set contained in R,{e)'' . If there exists 
to such that for every t < t' < to, Ct C Cf and lime(C) = B, then Ext(i?, i?(e)) 
has the same homotopy type as C . 

Proof. Hardt's Triviality theorem implies that there exists to > Oj and a homeo- 
morphism 

4>to '■ Cto ^ (07*0] ^ Uo<t<(oCt 

which preserves Cto ■ Replacing to by e gives a homeomorphism 

0e : C X (0,e] ^ Uo<t<eCt. 

Defining 

by 

V'(a;, s) = 7ri...fe o s), if s > 
V'(a;,0) = lims^o+ T^i...k o0(.t, s), 
it is clear that ip is a semi-algebraic retraction of C to Ext(i?, R(£}). □ 

We now prove Proposition 4. 

Proof, (of Proposition 4) Apply Lemma 5.1 to Co- and Ext(7o(7?.(cr-)), R2j-i): thus 
Co can be semi-algebraically retracted to Ext(7CT(7?.(cr_)), R2j-i). 

Since Ext(7o(7?.((T_)), R2j-i) is semi-algebraically contractible, so is Co. □ 

We now prove that the sets Ext (Co, R') form a covering of Ext (5*, R'). 
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Proposition 5. (Covering property) 

Ext(S',R')= U Ext(C^,R'). 

o-GSign(^,S) 

The proposition is an immediate consequence of the following stronger result. 
Proposition 6. 

Ext(S',R')= IJ 7^(CT^,R"^■). 

<TGSign(>t,S) 

Proof. By definition, Ext(5, R') D |J n{a+,R"'). We now prove the reverse 

aGSign{A,S) 

inclusion. Clearly, wc have that S = (J Ti{a,R). Let x G Ext(S', R') and a 

aeSign{A,S) 

be the sign condition of the family ^ at a; and let level((T) ~ j. If x G 7^((t1, R''^), we 
are done. Otherwise, there exists P & A, such that x satisfies either < P{x) < e2j 
or -£2j < P{x) < 0. Let e = {P G ^ | lim^^^- P{x) = 0}. Clearly #B = f > j. 
Let y ~ limj2j x. Since, Ext(S', R') is closed and bounded and x G Ext(S', R'), y is 
also in Ext(S', R'). Let t be the sign condition of ^ at y with level(r) = j' > j. 
If X G 7?,(t^,R' ) we arc done. Otherwise, for every P £ A such that P{y) = 
0, we have that —£2/-! < P{x) < e-zj'-i, since limg^j (P(a;)) = P{y) ~ and 
£2j'-i ^ £2j- So there exists P £ A such that x satisfies either < P{x) < e2j' 
or — £2j' < P{x) < 0, and we replace B hy {P £ A \ hnig^^, P{x) = 0}, and y by 
y = linie^^, X. This process must terminate after at most t steps. □ 

Algorithm 2. (Covering by Contractible Sets) 
Input. 

— a finite set of s polynomials V C D[Xi, . . . ,Xk\ in strong fc- general 
position on R'^, with deg(i-i) < d for 1 < i < s, 

— a T'-closed semi-algebraic set S, contained in the sphere of center 
and radius r, defined by a 'P-closed formula (j). 

Output, a set of formulas {(pi, . . . , (f>M} such that 

— each 7^(0i,R''") is semi- algebraically contractible, and 

— y 7e(0j,R"') = Ext(5',R'). 

l<z<M 

Procedure. 

Step 1 Let Q = Xf + . . . + Xl-r'^. Call Algorithm 1 (Parametrized Bounded 
Connecting) with input Q^V. Let A be the family of polynomials 
output. 

Step 2 Compute the set of realizable sign conditions Sign(^, S) using Algo- 
rithm 13.37 (Sampling on an Algebraic Set) in [6]. 

Step 3 Using Algorithm 14.21 (Quantifier Elimination) in [6], eliminate one 
variable to compute the image of the semi-algebraic map 7cr_ . Finally, 
output the set of formulas {0c, | a G Sign(^, S*)} describing the 
semi- algebraic set Ca- 
Complexity. The complexity of the algorithm is bounded by 5*'^+^^ (P^^''\ 
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Proof of correctness. The correctness of the algorithm is a consequence of 
Propositions 4 and 5 and the correctness of Algorithm 1 (Parametrized Bounded 
Connecting), as well as the correctness of Algorithms 13.37 and 14.20 in [6]. □ 

Complexity analysis. The complexity of Step 1 of the algorithm is bounded 
by s'^'^^dP^'^ \ where s is a bound on the number of elements of P and c? is a 
bound on the degrees of the elements of V, using the complexity analysis of Al- 
gorithm 1 (Parametrized Bounded Connecting). The number of polynomials in A 
is s'^'^^cf^'^ ' and their degrees are bounded by d'-"-'^ \ Thus the complexity of 
computing Sign(^, S) is bounded by s('=+i) d'^^'' ^ using Algorithm 13.37 (Sam- 
pling on an Algebraic Set) in [6]. In Step 3 of the algorithm there is a call to 
Algorithm 14.21 (Quantifier Elimination) in [6]. There are two blocks of variables 
of size k and 2 respectively. The number and degrees of the input polynomials 
are bounded by s'^+^d'^^*-' ) and d'-''^'^ ' respectively. Moreover, observe that even 
though we introduced 2s infinitesimals, each arithmetic operation is performed 
in the ring D adjoined with at most 0{k) infinitesimals since the polynomials 
{P, P ± e2j, P ± £2^-1 J P&PA<j^s} are in strong general position. Thus, the 
complexity of this step is bounded by 5^'^+^' d'-"-'' ^ using the complexity analysis 
of Algorithm 14.21 (Quantifier Elimination) in [6] and the fact that each arithmetic 
operation costs at most d'^''^ ' in terms of arithmetic operations in the ring D. □ 

6. Topological Preliminaries 

Wc first recall some results from algebraic topology which enable us to compute 
the first Betti number of any given closed and bounded semi-algebraic set, from the 
inclusion relationships amongst the connected components of the various pair-wise 
and triple- wise intersections of the elements of a covering of the given set by a finite 
number of contractiblc sets. 

6.1. Generalized Mayer- Vietoris exact sequence. Let Ai,...,An be sub- 
complexes of a finite simplicial complex A such that A = Ai U ■ ■ ■ U An- Note that 
the intersections of any number of the sub-complexes, A^, is again a sub-complex 
of A. We will denote by Aig,,,,Sp the sub-complex Aig n ■ • ■ n Ai^. 

Let C*(A) denote the Q- vector space of i co-chains of A, and C*{A), the complex 

y C-^A) C^A) C*+i(A) ^ • • ■ 

where d : C''{A) ^^^{A) are the usual co-boundary homomorphisms. More 
precisely, given to € C'^{A), and a q + 1 simplex [ao, . . . , Og+i] G A, 

(6.1) duj{[ao,. . . ,aq+i]) = ^ (-l)'w([ao, . . . , a^, . . . , a^+i]) 

0<i<ij+l 

(here and everywhere else in the paper " denotes omission) . Now extend duj to a 
linear form on all of Cq+i{A) by linearity, to obtain an element of C"^+^(A). 
The generalized Mayer- Vietoris sequence is the following: 

^ c'iA) ^ n^'(^''") n c'iA,„,j 

io io<il 

■■■'-^ n n c-(A.„,...,^,j... 

io<-<ip io<---<ip+i 
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where r is induced by restriction and the connecting homomorphisms 5 are described 
below. 

Given an w S njo< - <i ^'(^'o ip) define 5{ijj) as follows: 

first note that 5{lu) € Y\i„<---<i^^i C'(^io,...,ip+i): ^-^id it suffices to define (S(a;)io,...,ip^ 
for each (p + 2)-tuplc < zq < ■ • • < ip+i < n. Note that, 5{ijj)iQ^,,,^i^_^-^ is a linear 
form on the vector space, Cq{AiQ,,,,^i^^^), and hence is determined by its values 
on the g-simplices in the complex ^^(,,...,4^4.1. Furthermore, each g-simplex, s S 
^io,...,ip+i is automatically a simplex of the complexes A^^ ^- ^ , < i < p + 1. 

We define, 

(M.„,....p+.(^)= E (-i)H,...,A,....p..(^)' 

0<i<p+l 

(here and everywhere else in the paper " denotes omission). The fact that the 
generalized Mayer- Vietoris sequence is exact is classical (see [2] for example). 

The cohomology groups H'^{AiQ^,,,^i^) are isomorphic to the Q- vector space of 
locally constant functions on and the induced homomorphisms, 

5^:H*{A.,,_,^)^H*{A,,_,^^,) 

arc then given by generalized restrictions, i.e. for 

<^ S ©l<io<---<ip<n^°(^io,...,ip)i 

a locally constant function on 

p 

Sp{<|>)^o,...,^,^, = E("l)'<,....»;,....*P+J^-o......p+i- 

1=0 

The following proposition provides the key tool for computing the first Betti 
number. 

Proposition 7. Let Ai, . . . ,An be sub-complexes of a finite simplicial complex A 
such that A = ^1 U • ■ • U An and each is acyclic, that is H'^{Ai) = Q and 
Hi{A,) = for all g > 0. Then, 61(A) = dim(Ker((52)) - dim(Im((5i)), with 

i i<j i<j<£ 

To prove Proposition 7, we consider the following bi-graded double complex 
with total differential D = 5+ (-l)Prf, where 



io,...,ip 
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d 

V 
V 

V 





d 




d 




There are two spectral sequences (corresponding to taking horizontal or vertical 
filtrations respectively) associated with AiP''^ both converging to H'^{M.). The first 
terms of these are 'Ei = HsM,'E2 = HdHsM, and "Ei = HdM,"E2 = HsHdM. 
Because of the exactness of the generalized Mayer- Vietoris sequence, we have that, 



and 



'E,= 












d 

]" 

C\A) 
d 

C°(A) 









jo jo to jo 





'En 



H^{A) 

H^{A) 

H^{A) 

H°{A) 



The degeneration of this sequence at E2 shows that H'^[M.) = H*{A). 
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The initial term " Ei of the second spectral sequence is given by, 





s ^ 


nio<ii 










s 


T\io<ii 






.J 




s ^ 


nio<ii 


H\A,„,,,) 








s ^ 


nio<ii 




W . . H°(Ai„ „ 





The cohomology groups if°(Ai(,,...,ip) occurring as summands in the bottom row 
oi" El are isomorphic to the Q-vector space of locally constant functions on A^^ i^ 
and the homomorphisms, "di : — » "j^p+^'O are then given by generalized 
restrictions, i.e. for 



S ©l<io<---<ip<ri^f°(^io,...,ip)i 



with each 



K->,---,ip+i ^ locally constant function on Ai^^ . ^i 
p 



i,...,ip+i l^io- 



40,---,«i,---,4p + 



Proof. (Proof of Proposition 7) Since, H'^{Ai) = for all q > 0, all the terms in the 



first column of "Ei are zero except the bottom term, and clearly "(ij'^ = 0. Thus 



E'^'' and "K'^' = 0. Thus, H\A) ^ "E, 



□ 



7. Computing the first Betti number in the "P-closed case 



Let S" be a "P-closed semi- algebraic set. We first replace 5" by a 'P*-closed and 
bounded semi-algebraic set, where the elements of V* are slight modifications of the 
elements of V, and the family V* is in general position and hi{S*) ~ bi{S),0 < i < k. 

Define 

i<i<fc 

where d! is the smallest number strictly bigger than the degree of all the polynomials 
in V . Using arguments similar to the proof of Proposition 13.7 in [6], it is easy to 
see that the family V* of polynomials Pi — 5Hi, Pi + SHi, with Pi G P. is in general 
position in R((5)'^. 

Lemma 7.1. Denote by S* the set obtained by replacing any Pi > in the definition 
of S by Pi > —SHi and every Pi < in the definition of S by Pi < SHi. If S is 
bounded, the set Ext(S', R((5)*' is semi- algebraically homotopy equivalent to S* . 

Proof. Note that S is closed and bounded, lim^ 5* = S',and St C Sf. The claim 
follows by Lemma 5.1. □ 



Algorithm 3. (First Betti Number of a 'P-closed Semi-algcbraic Set) 



Input. 
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— a finite set of poiynomials V C D[Xi, . . . , Xk], 

— a formula defining a T'-closed semi-algebraic set, S. 
Output. 61(5*). 

Procedure. 

Step 1 Let £ be an infinitesimal. Replace S by the semi-algebraic set T defined as 
the intersection of the cylinder S x R(e) with the upper hemisphere defined 
by e\Xl + ... + XI + Xl^,) = 1, Xfe+i > 0. 

Step 2 Replace T by T* using the notation of Lemma 7.1. 

Step 3 Use Algorithm 2 (Covering by Contractible Sets) with input e'^{X^ -I- . . . -|- 
Xl + — 4 and P*, to compute a covering of T* by closed, bounded 

and contractible sets, Ti, described by formulas (jii. 

Step 4 Use Algorithm 16.27 (General Roadmap) in [6] to compute exactly one 
sample point of each connected component of the pairwise and triplewise 
intersections of the T^'s. For every pair i,j and every k compute the inci- 
dence relation between the connected components of Tj*^, and T*j as follows: 
compute a roadmap of T^*, containing the sample points of the connected 
components of T*-^ using Algorithm 16.27 (General Roadmap). 

Step 5 Using linear algebra compute bi{T*) = dim(Ker((S2)) — dim(Im((5i)), with 

Complexity. The complexity of the algorithm is bounded by (sd)*^"'^' , 
where s = and d = maxpgp deg(P). 



Proof of correctness. First note that T is closed and bounded and has the 
same Betti numbers as 5*, using the local conical structure at infinity. It follows 
from Lemma 7.1 that T and T* have the same Betti numbers. The correctness 
of the algorithm is a consequence of the correctness of Algorithm 2 (Covering by 
Contractible Sets), Algorithm 16.27 (General Roadmap) in [6], and Proposition 7. 
□ 

Complexity analysis. The complexity of Step 3 of the algorithm is bounded 
by s(fe+i)'dO(fe«) 

using the complexity analysis of Algorithm 2 (Covering by Con- 
tractible Sets) and noticing that each arithmetic operation takes place a ring con- 
sisting of D adjoined with at most k infinitesimals. Finally, the complexity of Step 
4 is also bounded by (sd)*^"'^' , using the complexity analysis of Algorithm 16.27 
(General Roadmap) in [6]. □ 

8. Replacement by closed sets without changing homology 

In this section, we describe a construction due to Gabrielov and Vorobjov [12] 
for replacing any given semi-algebraic subset of a bounded semi-algebraic set by 
a closed bounded semi-algebraic subset and strengthen the result in [12] to prove 
that the new set has the same homotopy type as the original one. Moreover, the 
polynomials defining the bounded closed semi-algebraic subset are closely related 
(by infinitesimal perturbations) to the polynomials defining the original subset. In 
particular, their degrees do not increase, while the number of polynomials used in 
the definition of the new set is at most twice the square of the number used in the 
definition of the original set. 
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Let C C R[Xi, . . . ,Xk] be a finite set of polynomials with t elements, and let 
S" be a bounded C-closed set. We denote by Sign(C, S) the set of reafizable sign 
conditions of C whose realizations are contained in S. 

Recall that, for a € Sign(C) wc define the level of a as #{P e C\(7{P) = 0}. As 
before let, e2t 3> £2t~i ^ ■ • • ^ £2 3> > be infinitesimals, and wc will denote 
by Rj the field R(e2t) • ■ • {£%) and denote by R' the field Ri. For i > 2t, R^ = R and 
for i < 0,Ri = R'. 

We now describe the construction due to Gabrielov and Vorobjov. For each level 
m,0 < m < t, we denote by Sign,„(C,S') the subset of Sign(C,S') of elements of 
level m. 

Given a £ Sign,„ (CjS") denote by TZ{a'^) the intersection of Ext(S', R2m) with 
the closed semi-algebraic set defined by the conjunction of the inequalities, 

~£2m < -P < £2m for cach P £C such that <j{P) = 0, 
P > 0, for each P e C such that cr(P) = 1, 
P < 0, for each P e C such that cr(P) = -1. 

and denote by, TZ{ij°i_) the intersection of Ext(S', R2m-i) with the open semi-algebraic 
set defined by the conjunction of the inequalities, 

— £2m-i < P < e2m-i for cacli P e C such that cr(P) = 0, 
< P > 0, for each P e C such that a{P) = 1, 
P < 0, for each P e C such that a{P) = -1. 

Notice that, 

Ext(7^(a),R2™) C7^(a;), 
Ext(7^((T),R2™_l) C7^«). 

Let X <Z Shea. C-semi-algebraic set such that X = (J ^(c) with S C Sign(C, S). 

creS 

We denote = S n Sign„j(C,S') and define a sequence of sets, X™ C R"^, 
< m < i inductively. 

• Define X° = Ext(X,R'). 

• For < m < t, wc define 

X"+i=[x'"U IJ Ext(7Z(a;),R') I \ |J Ext(7^(cr° ), R') 

V CTeS„ / creSign„^(C,S)\S„ 

We denote by X' the set 

The following theorem is a slight strengthening of a result in [12] (where it is 
shown that the sum of the Betti numbers of X and X' are equal) and the proof is 
very close in spirit to the one in [12]. 

Theorem 8.1. The seis Ext (X, P') and X' are semi- algebraically homotopy equiv- 
alent. In particular, 

H,{X)'^H,[X'). 

For the purpose of the proof we introduce several new families of sets defined 
inductively. 

For each p, 0<p<t + lwc define sets, Yp c rL, Zp c RL-i as follows. 
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• We define 

y/ = Ext(X,R2p)U U 7^(a;), 

zp = Ext(y/,R2p_i)\ U 7^«). 

f76Signp(C,S)\Sp 

• For p < m < t, we define 

y;"+i= [y;»u y Ext(7^«),R2p)) \ (J Ext(7^(a°),R2p) 

\ o-eS„ / creSign„(C,S)\E,„ 

^p"+'=[^™U y Ext(7e(a;),R2p-i)) \ y Ext(7e(a;),R2p-i). 

We denote by C R2p (respectively, Zp C R2p-i) the set Y^^^ (respectively, 

Note that 

• X ^ Yt+i = Zt+i, and 
. Zo = X'. 

Notice also that for each p,0 < p < t, 

(1) Ext(ZpP+i,R2p)cr/, 

(2) ZP cExt(r/',R2p_i). 

The following inclusions now follow directly from the definitions of Yp and Zp. 

Lemma 8.2. For each p,0 < p < t, 

(1) Ext(Zp+i,R2p) C Yp, 

(2) Zp cExt(yp,R2p-i). 

We now prove that in both the inclusions in Lemma 8.2 above, the pairs of sets 
arc in fact semi- algebraically homotopy equivalent. These suffice to prove Theorem 
8.L 

Lemma 8.3. For 1 < p < t, Yp is semi- algebraically homotopy equivalent to 
Ext(Zp+i, R2p). 

Proof. Let Yp{u) C R2p+i denote set obtained by replacing the infinitesimal e2p in 
the definition of Yp by u, and for uq > 0, we will denote by Let i^((0, uq]) C R2^+i 
the set {{x,u)\x € Yp{u),u € (0,mo]}. 

By Hardt's triviality theorem there exist uq G R2p+i, uq > and a homeomor- 
phism, 

ip : Yp{ua) X (0, uq] Yp{{0, uo]), 

such that 

(1) 7rfc+i((/)(x,u)) = u, 

(2) -0(a;,uo) ~ (x^uq) for x £ Yp{uo), and 

(3) for all u G (0, uq], and for every sign condition a of the family, Up^clP, -Pi 
£2*, . . . , P" ± e2p+i}, ''Pi', u) restricts to a homeomorphism of 7?.(cr, yp(uo)) 
to n{a,Yp{u)). 



20 



SAUGATA BASU, RICHARD POLLACK, AND MARIE-FRANgOISE ROY 



Now, we specialize uq to e2p and denote the map corresponding to ip by (j)- 
For cr g Sp, we define, TZ{cr'^_^) to be the set defined by, 

{-2e2p < P < 2e2p, for each P eC such that cr(P) = 0, 
P > —e2p, for cacli P £ C sucli that cr(P) = 1, 
P < £2p, for each P & C sucli that cr(P) = —1. 

Let A : ^ R2p be a semi- algebraic continuous function such that, 

A(x) = 1, on Fp n U^eE^7^(a;), 
A(x) = 0, on Fp \ U,eE^7^(a;+), 
< A(.t) < 1, else. 

We now construct a semi-algebraic homotopy, 

h:YpX [0,e2p] Yp, 

by defining, 

h{x,t) = 7ri...fc o A(a;)< + (1 - A(a;))e2p), 
for < t < £2p, 
0) = liuit^Q^ h{x , t) , else. 

Note that the last limit exists since S is closed and bounded. We now show that, 
h{x,0) G Ext(Zp+i,R2p) for all x eYp. 

Let X GYp and y = h{x, 0). 

There are two cases to consider. 

A(x) < 1: In this case, x G Ext(Zp-|_i, R2p) and by property (3) of ((> and the fact that 

A(.t) < 1, y e Ext(Zp+i,R2p). 
A(a;) = 1: Let ay be the sign condition of C at y and suppose that y ^ Ext(Zp+i, R2p). 
There are two cases to consider. 
Gy e S: In this case, y & X and hence there must exist r G Sign„j(C, S) \ !],„, 

with m > p such that y £ TZ{t°). 
Uy ^Yj-. In this case, taking t = Uy, level(T) > p and y £ TZ{t!^). 

It follows from the definition of y, and property (3) of 0, that for any m > p, 
and every p G Sign„(C, 5*), 

— 2/ G 7^(/0+) imphcs that x G 7?.(/9° ), and 

— X G Tl{p%) implies that y G TZ{p'^). 
Thus, X ^ Yp which is a contradiction. 

It follows that, 

(1) h{-,S2p) : — > is the identity map, 

(2) hiYp,0) = Ext(Zp+i,R2p), and 

(3) h{-, t) restricted to Ext(Zp+i, R2p) gives a semi-algebraic homotopy between 
/^(■,e2p)lExt(Zp+i,R2p) = i^Ext(Zp+i,R2p) ^'^'^ ^(■'0)lExt(Zp+i,R2p)- 

Thus, Yp is semi-algebraically homotopy equivalent to Ext(Zp-|_i, R2p). □ 

Lemma 8.4. For each p, < p < t, Zp is semi-algebraically homotopy equivalent 
to Ext(yp,R2p-i). 

Proof. For the purpose of the proof we define the following new sets for u G R2p- 
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(1) Let Zp{u) C be the set obtained by replacing in the definition of Zp, 
62 j by 62 j — u and £2^-1 by £2^-1 + u for all j > p, and £2p by e2p — u, and 
£2p-i by u. For uq > we will denote by Zp((0,uo]) the set {{x^u) \ x £ 
Zp{u),u e (0,uo]}. 

(2) Let Y.p{u) C be the set obtained by replacing in the definition of Yp, 
£2j by £2j — u and £2^-1 by £2j~i + u for all j > p and S2p by by £2p — u. 

(3) For a e SIGN™(C,S'), with m > p, let 7Z(cr^)(u) C denote the set 
obtained by replacing e2m by £2m — u in the definition of TZ{a'^). 

(4) For cr e SIGN„i(C,5'), with m > p, let n{a^){u) C Rap denote the set 
obtained by replacing £2m-i by £2m-i + u in the definition of 7^((t° ). 

(5) Finally, for a € SIGNp(C, S) let 7^(cr^)(^i) C R2P-1 denote the set obtained 
by replacing in the definition of TZ{(7^), £2p-i by u. 

Notice that by definition, for any u,v € R2p with < u < v, Z'p{u) C Yp{u), 
Z;{v) C Z;,{u), Y;{v) C F;H, and 

U r;(s)= U z;{s). 

0<s<u 0<s<u 

We denote by Z^ (respectively, Yp) the set Zp(£2p_i) (respectively, l^'(£2p-i)). It 
is easy to see that Y^ is semi- algebraically homotopy equivalent to Ext (Yp, R2p_i), 
and Zp is semi- algebraically homotopy equivalent to Zp. We now prove that, Y^ is 
semi- algebraically homotopy equivalent to Zp, which suffices to prove the lemma. 

Let ji : Yp —> R2p-i be the semi-algebraic map defined by 

fi{x) = sup {u \ X e Zp{u)}. 

M6(0,e2p-i] 

We prove separately (Lemma 8.5 below) that fi is continuous. Note that the 
definition of the set Zp{u) (as well as the set Yp{u)) is more complicated than the 
more natural one consisting of just replacing £2p-i in the definition of Zp by u, 
is due to the fact that with the latter definition the map fj. defined below is not 
necessarily continuous. 

We now construct a continuous semi-algebraic map, 

/i:yp'x[o,£2p-i]^rp' 

as follows. 

By Hardt's triviality theorem there exist uq G R2p, with uq > and a semi- 
algebraic homeomorphism, 

i^:Z;(uo)x(0,uo]->^;((0,«o]), 

such that 

(1) TTk+lii^ix^u)) = U, 

(2) 1/1(2;, uo) = {x,uq) for X G Zp{uo), and 

(3) ip{-,u) restricts to a homeomorphism of TZ{a, Zp{uo)) to TZ{a, Zp{u)), for 
every sign condition a of the family, Up^c{P, ± £2*, • • • , -P ± £2p+i}! for 
all u G (0, uq]. 

We now specialize uq to £2p-i and denote by (j) the corresponding map, 

0: X (0,£2p-l] ->Z;((0,£2p-l]). 
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Note, that for every u, < u < £2^-1, (f> gives a homeoniorphisni, (f>u ■ Z'p{u) Z'p. 
Hence, for every pair, w, it', < u < u' < £2^-1, we have a homeomorphism, 9u^u' ■ 
Zp{u) — > Zp(u') obtained by composing 0„ with For < u' < m < e2p-i, we 
let 9u u' be the identity map. It is clear that 0„ varies continuously with u and 
u'. 

For X G Yp, t G [0, £2p-i] we now define, 

h{x,t) = 6*^(2.), 

It is easy to verify from the definition of h and the properties of 4> listed above 
that, h is continuous and satisfies the following. 

(1) 0) : Fp' ^ Yp' is the identity map, 

(2) /i(r;,e2p-i) = ^;, and 

(3) h{-,t) restricts to a homeomorphism x t ^ for every t g [0,e2p-i]- 
This proves the required homotopy equivalence. □ 

We now prove that the function /i used in the proof above is continuous. 

Lemma 8.5. The semi- algebraic map fj. : R.2p-i defined by 

fj.{x) — sup {u \ X e Z'p{u)} 

Me(0,e2p-i] 

is continuous. 

Proof. Let < (5 <C £2p-i be a new infinitesimal. In order to prove the continuity 
of fi (which is a semi- algebraic function defined over R2p-i), it suffices to show that 

limExt(^,R2p-i((5))(a::') = limExt(^,R2p-i((5))(a;) 

5 S 

for every pair of points x,x' G Ext(yp', R2p-i ((5)) such that lim^ x = lim^ x' . 

Consider such a pair of points x,x' £ Ext(yp', R2p-i (5)). Let u £ (0,e2p-i] be 
such that X e Z'^iu). We show below that this implies x' € Zp{u') for some u' 
satisfying lim^ u' = lim^ u. 

Let m be the largest integer such that there exists a S with x G TZ{(t'^)(u). 
Since x G Zp(u) such an m must exist. 

We have two cases: 

(1) m > p: Let a G with x G TZ{a'^){u). Then, by the maximality of m, we 
have that for each P E C, cr{P) ^ implies that lim^ P{x) ^0. As a result, 
we have that x' G 7^(cr!j.)(u') for all u' < u — maxpgp o-(P)=o \P{x) — P{x')\, 
and hence we can choose u' such that x' G TZ{cr'^){u') and lim^ u' = lim^ u. 

(2) m<p: \ix' ^ Z'p(u) then since x' G Y^^ C Fp'(u), 

^' e U^eSignp(C,S)\Ep^«)(w)- 

Let (7 G Signp(C, S") \ Sp be such that x' G 7^((t° We prove by contra- 
diction that lim max = u. 

S PeV,a{P)=0 

Assume that 

lim max \P(x')\ ^ u. 

Since, x ^ TZ{a''^)(u) by assumption, and lim^ x' = lima x, there must exist 
P £ C, cr{P) ^ 0, and \misP{x) = 0. Letting r denote the sign condition 
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defined by t(P) = if \ims P{x) = and t{P) = <j{P) else, we liave tfiat 
level(T) > p and x belongs to both 7?.(t°)(u) as well as 7^(t^)(u). 

Now there arc two cases to consider depending on whether r is in E or 
not. If r S E, then the fact that x G TL{t'^){u) contradicts the choice of m, 
since m < p and lcvel(r) > p. If r ^ S then x gets removed at the level 
of T in the construction of Z'p{u), and hence x G TZ{p'^){u) for some p G S 
with level(p) > level(r) > p. This again contradicts the choice of m. Thus, 
lim max |P(a;')| = u and since x' ^ Uo-esiEn (c S)\s T^{<^+){'u') for all 

u' < max \P(x')\, we can choose u' such that lim^ m' = lim^ u, and 

x' ^ U^gsignp(c,s)\Sp7^(c^^;:)(^i'). 

In both cases we have that x' G Zp(u') for some u' satisfying lim^ u' = lim^ u, 
showing that lim5^(a;') > lim5^(a;). The reverse inequality follows by exchanging 
the roles of x and x' in the previous argument. Hence, lim^ /^(x') = \misp,{x), 
proving the continuity oi p,. □ 

Proof. (Proof of Theorem 8.1) The theorem follows immediately from Lemmas 8.3 
and 8.4. □ 



9. Computing the first Betti number of a general semi-algebraic set 

In this section wc describe the algorithm for computing the first Betti number 
of a general semi-algebraic set. We first replace the given set by a closed and 
bounded one, using the construction described in the previous section. We then 
apply Algorithm 3. 

Algorithm 4. (First Betti Number of a V- Semi-algebraic Set) 



Input. 

— a finite set of polynomials V C D[A"i, . . . , X^], 

— a formula defining a T'-semi-algebraie set, S. 
Output. &i(T). 

Procedure. 

Step 1 Let £ be an infinitesimal. Define S as the intersection of Ext(S', (e)) with 
the baU of center and radius l/er. Define Q as U {e^(^i -I- . . . -I- -I- 
^k+i) ~ 4: ^fc+i} Replace S by the Q- semi-algebraic set S defined as the 
intersection of the cylinder S x R(e) with the upper hemisphere defined by 
e\Xl + ... + XI + Xl^,) = 4, Xfe+i > 0. 

Step 2 Using the Gabrielov-Vorobjov construction described above, replace T by 
a Q'-closed set, T' . 

Step 3 Use Algorithm 3 to compute the first Betti number of T' . 

Complexity. The complexity of the algorithm is bounded by (sd)*^'''^', 
where s = #7-" and d = maxpgp deg(P). 



Proof of correctness. The correctness of the algorithm is a consequence of 
Theorem 8.1 and the correctness of Algorithm 3. □ 

Complexity analysis. In Step 2 of the algorithm the cardinality of Q' is 2(s + 1)^ 
and the degrees of the polynomials in Q! are still bounded by d. The complexity of 
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Step 3 of the algorithm is then bounded by (sd) 
of Algorithm 3. 



using the complexity analysis 



□ 



10. Computing Connected Components 



If one is interested in computing semi-algebraic descriptions of the connected 
components of a given semi-algebraic set, then using Algorithm 1 (Parametrized 
Bounded Connecting) it is possible to do so with a complexity making precise the 
one of previously known algorithms, whose complexities were of the form {sd)^ 
(see [16]). We have the following theorems (we refer the reader to [6] for details of 
the proof). 

Theorem 10.1. If7.{Q,Y&) is an algebraic set defined as the zero set of a poly- 
nomial Q € D [Xi , . . . , Xk] of degree < d, then there is an algorithm that outputs 
quantifier free formulas whose realizations are the semi- algebraically connected com- 
ponents of7j{Q,Y&). The complexity of the algorithm in the ring generated by the 
coefficients of Q is bounded by d^^^ ^ and the degrees of the polynomials that appear 
in the output are bounded by 0{d)^ . Moreover, if D = Z, and the bitsizes of the 
coefficients of the polynomials are bounded by t, then the bitsizes of the integers 
appearing in the intermediate computations and the output are bounded by rd^^'' ' . 

Theorem 10.2. LetV ^ {Pi, . . . ,PJ C D[Xi, . . .,Xk] withdcg{Pi) <d,l<i<s 
and a semi-algebraic set S defined by a V quantifier-free formula. There exists an 
algorithm that outputs quantifier-free semi- algebraic descriptions of all the semi- 
algebraically connected components of S . The complexity of the algorithm is bounded 
by s^^^d'-'^^ \ The degrees of the polynomials that appear in the output are bounded 
l,y (]P{k Moreover, if the input polynomials have integer coefficients whose bitsize 
is bounded by r the bitsize of coefficients output is d*^''"' 'r. 
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